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ABSTRACT 

Interferometers require accurate determination of the array configuration in order to produce 
rehable observations. A method is presented for finding the maximum-Hkehhood estimate 
of the telescope geometry, and of other instrumental parameters, astrometrically from the 
visibility timelines obtained from observations of celestial calibrator sources. The method 
copes systematically with complicated and unconventional antenna and array geometries, with 
electronic bandpasses that are different for each antenna radiometer, and with low signal-to- 
noise ratios for the calibrators. The technique automatically focusses on the geometry errors 
that are most significant for astronomical observation. We apply this method to observations 
made with the Very Small Array and constrain some 450 telescope parameters, such as the 
antenna positions, effective observing frequencies and correlator amplitudes and phase shifts; 
this requires only 1 h of CPU time on a typical workstation. 
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1 INTRODUCTION 

Accurate knowledge of observing frequency, projected baselines 
and other geometric parameters is essential if interferometer data 
are to be turned into useful images. The accuracy required is typ- 
ically a small fraction of the operating wavelength. For example, 
to keep phase errors below 10 degrees for an interferometer oper- 
ating at 30 GHz, it is necessary to know the baseline length to an 
accuracy of better than 0.3 mm. 

The standard way in which interferometer geometries are de- 
termined is along the following lines. First, direct measurements 
are made using, for example, tape measures, theodolites, and mea- 
surements of electrical lengths of electronic systems. Then, obser- 
vations of celestial calibrators are carried out. For example, the ob- 
servation by a particular pair of antennas of a bright calibrator of 
known RA and Dec gives the collimation phase of that pair (the 
difference in electrical length of the two arms of that interferome- 
ter). How this quantity varies, for example, with hour angle and for 
a series of calibrators at different declinations, is also investigated. 
These standardmethods^i^^ Moran & 

Swenson''(|^9^)7anTTi71hese'^ayrT'seri£rof corrections to the 
directly measured geometry is produced, with the aim of ensuring 
that the physical baselines (in wavelengths) used in subsequent data 
analysis are correct. 

This approach is of proven effectiveness for determining the 
geometries of interferometer arrays such as the Very Large Array 
(VLA). It is not adequate, howev er, for interfer o meter s such as the 
Very S mall Array (VSA, see e.g. [Watson et ak] ( |2002| ) 



Scott et al 



(2002)), which is designed for imaging the CMB on scales > 10 



arcmin at 30 GHz, and thus has baselines of length < 2 m. At first 
sight this may seem surprising, as VSA baselines (in both metres 
and wavelengths) are so small. From the point of view of geometry 
determination from astrometry, however, the VSA is diiferent from 
most other interferometer arrays in two respects. 

• The collecting area of VSA antennas is very small and so there 
are few suitable celestial calibrators. Only the brightest calibrators 
are of use and signal-to-noise is a major issue. 

• For most arrays like the VLA, the principal swivel axes of 
each antenna (e.g. altitude and azimuth, or polar and declination) 
intersect, at least nominally. If an axis is incorrectly aligned (e.g. 
the polar axis does not quite point to the Pole), then to first order 
this has no effect on the measured phase of a source, and its effect 
on source amplitude will manifest as a simple pointing correction 
which a standard 5-point observation (at expected source position, 
then offsets North, South, East then West) will provide. In the case 
of the VSA, the table elevation axis and the individual horn rota- 
tion axes (see Section ^ do not intersect by ~ 1 m. This design 
helps reduce important systematic errors and also reduces cost, but 
it makes astrometric geometry determination significantly harder. 
If an axis is misaligned, then in general a phase error is also intro- 
duced. The identification of such errors is exacerbated by the lack 
of signal-to-noise. 

There are also further difficulties associated with determining the 
geometry of the VSA, which can be present in other interferome- 
ters. 

• As a result of, for example, the compliance needed for fitting 
horn-reflectors to cryostats, the position of the electrical centre of 
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each horn will be slightly different, manifesting as a physical base- 
line correction and a pointing correction. 

• A large observing bandwidth makes it hard to achieve a flat 
band pass, so the effective radio frequency becomes antenna de- 
pendent. 

• Unconventional geometry makes it difficult to make a straight- 
forward interpretation of a set of geometry calibration observations. 

As a result of all these factors, phase (and amplitude) effects 
arise from complicated causes that are particularly hard to deter- 
mine because of the lack of very high signal-to-noise calibrators. 
We have therefore implemented a maximum-likelihood method to 
determine the geometry, and other telescope parameters, astromet- 
rically from the visibility timeline data obtained during a calibra- 
tion observation of a sufficiently bright unresolved radio source. In 
the geometric model, we include as variables all the important di- 
mensions and angles that require accurate determination, including 
the positions of the antennas, their effective observing frequencies 
and phase shifts on different baselines. In finding the best-fit ge- 
ometry, the maximum-likelihood method correctly allows for the 
signal-to-noise ratio of the observations (on the assumption that the 
instrumental noise is Gaussian); the procedure also automatically 
guards against attempting to determine unimportant geometry er- 
rors that do not in fact much affect astronomical observation, since 
these similarly will not much affect calibrator observations. 

In this paper, we focus on the method in its application to the 
VSA, although we point out the usefulness of the approach gener- 
ally to radio and optical interferometry. In Section |^ we present the 
general maximum-likelihood method for the geometry calibration. 
The VSA is introduced in Section ^ where we develop a model 
of the telescope as a function of the parameters to be calibrated. 
Using this model, the likelihood method is applied to real VSA ob- 
servations in Section ^. Finally, our conclusions are presented in 
Section ^ 



2 ASTROMETRIC GEOMETRY CALIBRATION 

Given a multi-dish interferometer consisting of Na separate an- 
tennas, with antenna positions denoted by dj (i = 1, . . . ,N„) rela- 
tive to some origin, one can measure the correlated signals of any 
two antennas along Nj, = jNa{Na — I) independent baselines Dj 
(j = 1, . . . ,A'i), given in some ground-based reference frame spe- 
cific to the telescope. The phase of the signal is modulated by the 
geometric path difference that the incident radiation experiences 
between the antennas. If the baseline is represented by a vector u 
in a coordinate system (Mv-coordinates) such that the component u 
of the vector u = {u,v,w) is measured in the east-west direction, 
V in the north-south direction and w in the direction of the anten- 
nas' pointing centre, the path difference is given by the projection 
of the direction cosines x — {x,y,z) of the source on the sky onto the 
baseline components, i.e. hy A = u x. Because the field of view is 
a section of a sphere, the third direction cosine is z = \/l—x^ ^y^- 
Assuming a small field size or a sufficiently compact source, z ~ I . 

The response of an interferometer to an extended source can 
be obtained by integration over the solid angle of the source. Let us 
denote the sky brightness distribution of the source at the observ- 
ing frequency v by I{x,y,v). One defines the complex (not fringe- 
rotated) visibili^_^^^^u^v^w^V^^^s^ Moran & 
Swenson ligg^l ) 

f{u,v,w,v) = l°° A{x,y,v)I{x,y,v)e^'""-^"-'+"y+''-^dxdy, (1) 



where A(x,j,v) is the directional antenna (power) response pattern 
or primary beam (normalised to unity at its peak). To simplify the 
notation, we denote the real part '^['U («)] of the visibility V («) by 
and similarly 3(«) for the imaginary part. 
For point sources at positions xi^^ from the pointing centre, 
whose brightness distributions are given by I(x) = 8{x — x^), 
the integral (|j|) for the complex visibilities V (m, v) reduces to 

3?(«,v) = £A(xi,v)SiCos('27l-«-Xi) , 

A- \ c / 

3(«,v) = Y.Mxk,v)Sk&m[2Tl^u-xi,y (2) 

k 

where the sum is over all relevant sources in the field of view and 
Xk and Sk denote, for the kth source, respectively the position from 
the pointing centre and the flux. 

As the Earth rotates, the a-vector corresponding to a base- 
line Dj varies with the time f; at which the visibility sample was 
measured. We will denote the wv-coordinates of the jth baseline at 
a time t, (or alternatively at hour angle Hi) by uij = u{Dj,ti), and 
the corresponding visibility by 1^^. The transformation R(8,/f,-,L) 
that relates a baseline vector Dj in the telescope frame to the uv- 
coordinates is a simple rotation given by 

Uij = R{5Mi,L)Dj, (3) 

which will generally depend on the geographic latitude L of the 
telescope and the declination 8 and hour angle Hi of the phase- 
tracking centre. 

For a given a model of the telescope, the signal timelines for 
an observation of a point source can be predicted and compared 
with the observed data. If the source is sufficiently bright, one may 
then use a maximum-likelihood approach to determine the antenna 
geometry and a large number of other telescope parameters. For 
Gaussian instrumental noise the maximum-likelihood solution for 
the parameters can be found by minimising the standard x^-misfit 
statistic between the predicted visibilities 1^^ and observed visibil- 
ities 'l^'^, namely 

XH-^U^^U) = Y"^'7 [(9?°(«,7)-9^'(«<7))' 
i-j 

+(30(«,v)-3''(«<7))']- (4) 

In this expression, the sums run over all A'f time samples and A'^ 
baselines. The weights on the observed visibilities are given 
by Wij = 1 /o?- assuming that the noises are independent, which is 
true for the thermal Johnson noises from each antenna. For a given 
observation, is a function of Np parameters {p/} (/ = 1 , . . . ,Np), 
including, for example, the antenna positions on the table {rf,}, the 
correlator amplitudes {Cj}, (constant) residual phase shifts {^j} 
and effective observing frequencies {v,}. By minimising the x^- 
function with respect to these parameters, it is possible to find the 
maximum-likelihood estimate of the telescope geometry. 

Finding the global minimum of the x^-function is numeri- 
cally quite challenging. The optimisation takes place in an Np- 
dimensional space and, as outlined for the VSA below, one typi- 
cally has Np ~ 450. Moreover, the X^ = constant hypersurfaces are 
rather uneven. In particular, there are many local minima, and a 
straightforward (e.g. gradient search) minimisation algorithm may 
never find the correct basin of attraction. The existence of local 
minima can be understood intuitively, since a minimum occurs 
whenever the fringe pattern is fitted correctly on only one of the 
baselines. Even though the false minimum can easily be recognised 
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from the poor fit on at least some of the timehnes, a straightfor- 
ward optimisation would stop there. Hence a successful fit requires 
a good initial estimate. Simulations show that for the VSA the ini- 
tial antenna positions have to be measured manually to an accuracy 
of at least a third of a wavelength, i.e. roughly 3 mm. This is feasi- 
ble with a steel tape measure. It should also be noted that there is a 
certain degeneracy in the problem. For example, a rescaling of the 
baselines is equivalent to a scaling of the observing frequency. 

There are a number of possible numerical minimisation algo- 
rithms to choose from. Some of them er nploy gradient inf ormation, 
such as conjugate gradient algorithms ( Press et al 1992 ). Deriva- 
tives of the x^-functional with respect to the telescope parame- 
t ers can be found in Ap pendix Q Simulated annealing methods 



(Kirkpatrick et al. 19831) try to find their way through local min- 



ima. In this application, however, the number of parameters is very 
large and the /^-surface does not have any significant large-scale 
slope, rendering a minimisation in a large number of dimensions 
very difficult even for the annealing method. Given a good ini- 
tial esti mate, we found a c omparatively simple algorithm, Powell's 
method (Press et al. 1992), to be very robust and reliable. 




ground 



East 



Figure 1. The VSA antennas are mounted on a table with a hinge in east- 
west direction, which tracks a source in elevation E. The individual horns are 
tilted at an angle T| with respect to the table and track in table 'azimuth' p. 



3 THE VERY SMALL ARRAY 

The VSA consists of an array of 14 horn-reflector antennas 
mounted at variable positions on a table that can be steered in el- 
evation. The antennas track individually about their horn axes. A 
metallic enclosure serves as a ground shield. The observing fre- 
quency can be varied between 26 and 36 GHz, with a receiver band- 
width of 1.5 GHz and an effective system temperature of 30-35 K. 
The telescope was designed to be able to observe in two different 
configurations using antennas of different sizes. 

(i) The compact array consists of conical horn reflector an- 
tennas each with an aperture diameter of 14.3 cm. The primary 
beam has a FWHM of about 4.6° at a frequency of 34 GHz. The 
CMB power spectrum is measured in 10 independent bins of width 
A£ ^ 80 from £^ 100- 900. The VSA observed in the compact 
configuration between August 2000 and August 2001. 

(ii) The extended array is a scaled version of the compact array. 
The aperture diameter is 32.2 cm, sampling multipoles £ ~ 350 — 
1900 with a resolution of A£ ~ 200. The VSA has been observing 
in the extended configuration since September 2001. 

The table on which the VSA antennas are mounted can be 
tilted about a hinge aligned along the east-west direction (Fig. |l|). 
The maximum elevation angle £ of the table is 70° and the mini- 
mum is 0°. The antennas have an inclination relative to the table 
of T| = 35° towards the south, allowing observations of up to 35° 
from the zenith to the north and south. Each antenna can individ- 
ually track a source east-west up to angles of (3 = 40° from the 
meridian. 

For the VSA, the rotation R from (^) is a function of the ta- 
ble elevation e. If the elevation has not been determined accurately, 
an error Ae is introduced. Misalignments of the table hinge can 
be taken into account by polar angles 6 and (|) describing the de- 
viation of the hinge from the east-west direction. Using the three- 
dimensional rotation matrices Rv((|)), Rv((|>) and Ry{<^) for right- 
handed rotations about the corresponding coordinate axes, R is 
given by 

R(8,/f,L,Ae,e,x) 

= R,(|-5)R,,(-//)R,(L-|)Ry(e)R,(x)Rv(e + Ae). (5) 



The actual arrangement of the antennas on the table for the 
compact configuration in September 2000 is shown in Fig. ^ (a). 
The resulting projected baselines during a 6-hour observation of a 
field at declination 5 = 22° are shown in panel (b). Since the tele- 
scope table is stationary with respect to the ground, during a track 
on a field the samples u lie on a series of curves (or t<v-tracks). 
The electronics of the receiving system of the VSA are illustrated 
in Fig. |3| An antenna receives a signal at the radio frequency Vrf- 
After mixing with the two local oscillators, which produce a com- 
bined frequency Vlo = Vlo,i +Vlo.2' the frequency of the signal 
is reduced to the (2nd) IF frequency Vip = Vrf — Vlo- I" order to 
measure both the real and imaginary parts of the complex visibility, 
the signal is split into sin- and cos-channels with a phase difference 
of 90 degrees, yielding 28 such channels for A'^ = 14 antennas. To 
avoid a loss of signal coherence due to changing path lengths dur- 
ing a tracking run, appropriate geometric path delays are inserted 
into each channel. The correlated signals from each pair of anten- 
nas are integrated to produce the real and imaginary parts of the 
Ni, visibilities. To minimise the number of required splits of each 
cos- or sin-channel, the real part of the visibility for a given baseline 
can be formed from either the cos-cos or the sin-sin coinbinations, 
whereas the imaginary parts may consist of cos-sin or sin-cos. A 
lookup table is used to determine which cos- and sin-channels cor- 
respond to real and imaginary parts of a given baseline. 

The attenuation and passbands of filters and amplifiers are not 
spectrally flat, resulting in slightly different effective IF frequen- 
cies ViF.cos and VjF sin for each antenna. Additional path delays in- 
troduced by the cables and electronics between antennas and cor- 
relators also affect the visibilities, and the amplitudes of the signal 
differ due to the different electronic gains of the correlators. There- 
fore the amplitudes cj'^' and cj'' of the real and imaginary parts 

and phase shifts and (|)^'' vary for each baseline j. However, 
tests show that these are effectively constant over timescales of ap- 
proximately one month on a given baseline. 

(r) 

The raw data have to be corrected for inserted path delays p ■ ^ 

and p^.'j, where / = 1,2 denotes the cos and sin-channels from dif- 
ferent antennas corresponding to the real and imaginary part of the 
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Figure 2. Top: the compact atrangement of the 14 VSA antennas on the 
table in September 2000. The rectangular area indicates the shape of the 
lower part of the VSA table, with the hinge on the lower end. The dashed 
lines above the horns show the extension of the antenna mount ( j^usholrn^ 



2001 ). Bottom: the in'-coverage of the VSA using the antenna arrangement 
shown in (a). The observed field is at a declination of 8 = 22°. The duration 
of the observation is 3 hours either side of transit. 



yth baseline. The path delays are introduced at the different IF fre- 
quencies v|p'^. 1 and vjp ^. J. Note that this notation with 2x Nj, dif- 
ferent symbols obscures the fact that there are just 2xNa different 
path delays and IF frequencies, as the path delays are inserted per 
antenna and not per baseline. 

Gathering together the effects described above, the predicted 
visibilities for the jth baseline for observations of a set of point 
sources at positions xj^ with fluxes 5^ are thus given by 



X cos 271 



(r) [r) 

(r) %J,2 (r) 



1st LO (24 GHz) 



'-Xjj- 

2nd LO (8 GHz) 




2nd antenna 



Figure 3. Simplified schematic of the electronics of the receiving system of 
the VSA for compact array observations at a centre frequency of 33 GHz. A 
signal coming from an antenna at the radio frequency Vrp is reduced to the 
(2nd) IF frequency Vip after mixing with the signals coming from the two 
local oscillators. It is split into sin- and cos-channels, and a geometric path 
delay is inserted into each channel separately. The correlated signals from 
pairs of antennas are integrated to produce the real and imaginary parts of 
the visibility. 
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Uj-Xk 



, (6) 



where v = Vip + Vlo. l + is the radio frequency corresponding 
to a given IF frequency Vip. The superscripts ^ indicate that the 
visibilities have been theoretically predicted. Given the calculable 
dependence of uj for each baseline on the time t, one can thus 
calculate the quantities Vfj appearing in the expression for y} given 



4 ANALYSIS OF VSA CALIBRATION OBSERVATIONS 

In order to constrain the geometry, an observation with a good 
signal-to-noise ratio at a sufficiently high elevation is required. 
There are three possible calibration sources for the VSA com- 
pact array that fulfil these requirements: Cas-A, Cyg-A and Tau-A. 
These are essentially unresolved point sources for the VSA com- 
pact array. In the plots presented below, we use a 4-h calibration 
observation of Cas-A as an illustrative example. The optimisation 
with respect to antenna positions, correlator amplitudes and phases 
and the IF frequencies was carried out using Powell's algorithm 
down to a relative convergence tolerance of lO^-' in the value of y}. 
The entire calculation required around I h of CPU time on a Sparc 
Ultra workstation. 
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Figure 4. Comparison of observed and predicted fringes on VSA baseline antenna 1/antenna 2, using an observation of Cas-A. The panels show the real and 
imaginary parts, the amplitude and the phase of the visibility respectively, as a function of hour angle (in hours). The theoretical fringe predicted from the 
optimised telescope parameters is plotted with a solid line and the observed fringe with a dashed line. The units on the y-axis are arbitrary. 



4.1 Timelines 



4.2 Tlie -surface 



As an illustration of the fit obtained, in Fig. q we plot the observed 
fringes for VSA baseline antenna 1/antenna 2. The four panels 
show respectively the real and imaginary parts of the visibility, the 
amplitude and the phase, each as a function of hour angle. Each 
panel shows the timelines for the observed fringes (dashed) and the 
fringes predicted assuming the best-fit telescope geometry (solid), 
obtained from the -minimisation. The scatter of the observed real 
and imaginary parts around the predicted lines indicates the noise 
level of the observation. 

The quasi-sinusoidal fringe pattern shows an increasing fringe 
frequency for the 1/2 baseline. Since the data have not been fringe- 
rotated, the phase is going rapidly t hrough cycles of 2t i. On this 
baseline, small quadrature errors (see Watson et al. ( 2002 )) and dif- 



ferences in the amplitudes on the two channels are seen as a slight 
variation of the visibility amplitude and as structure in the visibility 
phase. We note that the best fit is in excellent agreement with the 
observations and matches the observed amplitudes and phases very 
closely. 



Confirmation that the optimisation procedure has indeed fully con- 
verged can be seen in Fig. ^ |^ and ^. These plots show one- 
dimensional cuts through the /^-surface in parameter space along 
some of the parameter coordinate axes. In each case, the y-axis in- 
dicates the value of the x^-function, with an arbitrary normalisa- 
tion, as a function of the deviation of the parameter value from 
the best-fit value. All plots are centred around a local minimum 
in the /^-surface, indicating that the optimisation procedure has 
fully converged, and that the obtained reconstruction is indeed the 
maximum-likelihood estimate. There is, however, a strong caveat 
in interpreting these one-dimensional x^-contours: they are by no 
means representative of the accuracy with which the parameters 
can be determined, which would require a marginalisation over all 
other parameters rather than a simple cut. In fact, the position and 
shape of the one-dimensional contour can be changed significantly 
by simply varying one or a few of the other parameters. To ensure 
that a minimum has been found, all parameters have to converge 
simultaneously. 

Fig. M shows a cut along the x-position of VSA antenna 1. 
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Figure 5. The -surface as a function of the .v-position of VSA antenna 1. 
The units on the x-axis are metres from the optimised position, the v-units 
are arbitrary. Low values of correspond to the most likely solution. The 
minimisation has converged to the local minimum in the centre. A vertical 
line indicates the position of the manually measured value of the position. 
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Figure 6. The x~-surface as a function of the y-position of antenna 1. Com- 
pare Fig. ^ It is interesting to note that adjacent minima are nearly as pro- 
nounced as the central one. 
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Figure 7. The x^-surface as a function of the IF frequency on the cos- 
channel of antenna 1. The A-axis is in units of 12 GHz from the optimal 
frequency, i.e. 0.01 corresponds to 120 MHz. 



Figure 8. An image of Jupiter observed by the VSA. The geometry and 
other telescope parameters have been simultaneously calibrated by obser- 
vations of Tau-A and Cas-A. The dynamical range of the image is 500:1, 
indicating that the phase errors due to the calibration are indeed negligible. 



Minima in correspond to the most likely solution. As for all 
other plots, the minimum has been found. A vertical line indicates 
the original input value of the position, as was manually measured 
on the table. Fig. ^ shows a cut along the y-position for the same 
antenna. It is obvious from the figure that, in this case, the calibra- 
tion observation used is only just sufficient to distinguish between 
the true minimum and the adjacent minima either side, which are 
nearly as deep. As mentioned previously, to avoid being trapped in 
a false local minimum, one requires a sufficiently accurate initial 
estimate. Nevertheless, it is still possible to detect false minima, 
since the resulting positioning mistake of the order of 1 cm induces 
significant phase errors on at least one of the baselines of which the 
corresponding antenna is a part. 

As a further illustration. Fig. ^ shows a cut along the IF fre- 
quency of the cos-channel of antenna 1. Once more, we see that the 
position obtained by the Powell algorithm lies indeed at the mini- 
mum for this parameter. 



4.3 Imaging with the VSA 

The parameters found from the calibration observations can be ap- 
plied to the data reduction of subsequent observations. Fig. ^ shows 
a VSA map of Jupiter, which has been produced using a calibration 
that combined observations of both Cas-A and Tau-A. This high- 
lights that the geometry fitting procedure can easily exploit several 
data sets simultaneously by simply adding up the x^-functionals 
for all observations. The dynamical range of the image is an excel- 
lent 500:1, showing that the geometry calibration on other sources 
was sufficiently successft il to render ph a se err ors negligible. Other 
VSA maps are shown in Watson et al. (2002), and the primordial 
CMB maps obtained from VSA observations in the compact con- 
figuration are presented by Taylor et al, (2002). 
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5 DISCUSSION AND CONCLUSIONS 

In this paper, we have presented a method for astrometric geometry 
calibration of interferometers by performing a likelihood analysis 
of visibility timelines produced during an observation of a reason- 
ably bright radio point source. The method has been demonstrated 
for the specific case of the Very Small Array, for which one must 
simultaneously constrain some 450 telescope parameters, such as 
the antenna positions, the effective observing frequency of each 
antenna which can vary due to different bandpass characteristics, 
and phase shifts on different baselines introduced in the electronic 
system or by different cable lengths. In particular, analysis of VS A 
observations of Tau-A and Cas-A sources show that the relative 
positions of the antennas on the telescope table can be determined 
to a small fraction of a wavelength. The quality of the fit is suffi- 
ciently good to enable high dynamic-range mapping with the VSA. 
This astrometric calibration is very straightforward and cheap, as 
it requires no equipment other than the telescope itself. Multiple 
calibration observations can easily be combined to produce a better 
parameter estimate by summing the x^-functionals for all available 
data sets. 

The calibration method is currently in active use by the VSA 
project and has successfully determined the geometry for several 
different antenna configuration since the start of its observing pro- 
gramme. The calibration of the VSA is described in more detail in 



Watson et al. 



(2002), and power spectra derived from the first year 



of ob^servationsjii^he^om^^ 
et al. ( |2002[ ). 

Although the astrometric calibration method has been applied 
specifically to the VSA, it can easily be adapted to other telescopes. 
It is worth noting, however, that there is one feature of the VSA 
which does not necessarily generalise to other CMB telescopes, 
even though it is common to most ground-based arrays. Each VSA 
antenna tracks the source separately in 'azimuth' {semi-comounted 
design). This has the advantage of providing a valuable check on 
systematic contaminations, because the path delays on a particu- 
lar baseline change during the observation and the sky signal can 
be recognised by its char acteristic fringe fr equency. On the Cos- 
mic Background Imager ( Padin et al. II2OO2I) and the Degree An- 
gular Scale Interferometer (Leitch et al. 2002), the antennas are 
mounted on a rigid tracking platform supported by an alt-azimuth 
mount. This design keeps the baseline projections with respect to 
the source constant and reduces the fringe rate on all baselines to 
zero, making the described geometry fitting impossible. A similar 
method can still be applied, however, if the fringe rate is artificially 
changed. For instance, the telescope could be kept at a fixed point- 
ing position while the source is within the antenna beam, thus pro- 
viding a non-zero fringe rate. 
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APPENDIX A: DERIVATIVES USED IN THE VSA GEOMETRY CALIBRATION 

For reference purposes, in this Appendix we quote the expressions for the derivatives of the x^-functional (^) combined with the visibilities 
used in the maximum-likelihood algorithm for the the VSA geometry calibration. 

The gradient is required by some numerical minimisation algorithms that use gradient information to accelerate the computation. Al- 
though not required by the minimisation algorithm, the curvature can be useful for some iterative minimisation algorithms and for quantifying 
errors or uncertainties on the reconstruction. Assuming that the maximised estimate of the parameters is indeed in the correct basin of attrac- 
tion on the likelihood surface and that this surface can be reasonably approximated by a Gaussian around the maximum, the uncertainties in 
the reconstruction can be obtained from the inverse of the Hessian matrix of j} . 



Al First derivatives of j} 

From (W), the x^-gradient with respect to the telescope and correlator parameters pi is 



9k 



N,.Ni, 

21' 



9lP(«,v)-9lO(«,v: 



d9?P(« 
Ap, 



'^' + (3-(«,7)-3"(«,V, 



d3P(«,- 
dp, 



where again the symbol 9l^(«,y) denotes the real part of the visibility predicted for the /th time sample on the j\h baseline. The dependence 
of the visibilities on the parameters can be explicit, as in the case of the correlator amplitudes, or implicit in the position of the samples «y , 
as for the antenna positions: 

d9^P(«,v) _ 39^P(«,7) ^ m^jujj) dujj ^^^^ 
dpi dp, duij dp, ■ 



■('•) 



We now substitute the expressions (|6|) appropriate for the VSA. For the sake of notational brevity, we define the arguments 
7^ (J'-) vM. ,1 



-Uj -Xk 



(r) (r) 1 

■^IF.I (r) %,2 (r) I , .(r) 

—p\ -—Pi } +r' 



(A2) 



and similarly for and v^'' . Then the derivatives of the predicted visibilities with respect to the Mv-loci, 
99lP(«) 



9u 

^ = 2;rC«£A(x,,vW)5,cos(2:rV'>)^x, 



('•) 



-2nd''^ £A(xi,v(''))Si sin ( 27t\|/(''' ) —Xk 



(A3) 



only depend on table parameters, whereas 
^ = I:A(.„v«)5,cos(2.V'- 
93P(«) 

d3i^{u) 

d3''(«) 



= £A(xi,v«)5isin 27tV' 



-2i:&^ £A(xi,v(''')5t sin (2 V' 



27lC('' £A(xi,v('))5iCos C2n\|/('' 
k ^ 



(A4) 



(A5) 



only depend on parameters of the electronic system. 
The derivatives with respect to the frequencies are 



99lP(«) 



9v 



-2iid'"> Y,MxkM''^)Sk sin hn^^i' 



IF. I 



99lP(«) 

w 

IF,2 



av 



-2nC(''' Y,A{xkM''^)Sk sin (2 V' 
k 



2c' 



(r) ('•) 

V2 Pi 



Uj-Xk, 



^^\uj-X,, 



(A6) 



and analogously for the imaginary parts. These derivatives apply to a single baseline only. In practice, a given frequency channel is used in 
several baselines. In order to compute the x^-derivatives one thus has to sum over all baselines containing the respective frequency channel. 
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The derivatives of the in'-positions with respect to the geometric parameters are given by 

3 



Lm^Di ifp,e{A£,e,0} 



otherwise 

where subscripts denote the vector components, Z) is a baseline vector and R(5,//;,L,Ae,9,x) is the transformation matrix from table to 
liv-coordinates introduced in (^): 

Uij = R{8,H,L,Az,Q,x)Dj({di}). 

Derivatives with respect to parameters other than the table parameters will vanish since they do not affect the geometry. Furthermore, a given 
baseline depends only on the positions of those two horns of which it is composed, and the resulting matrix is considerably sparse. 



A2 Curvature of 

Denoting the real part of the visibility predicted at the ith time sample on baseline j by the x^-curvature is given by 



2^2 



dpidpn 



'J 



?>Pn 



9m 



9t 



dpidpn 



+ 



du dp„ 
■^i j ^ •■'i j 



dpi 9h 3/7/ 



9pn 



33)} 9„/A (d% 33)} duij 



du dpi 
( 9« 



+ 



V 3p/i 



[v„V„3fy 



3"/7 

dpi 



(duij\ 



d'ui 



^ 

du dpidpn 



d% d^u. 



du dpidpn 



Again the expressions (A3)-(A6l can be used to substitute for the partial derivatives. 



